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Qs I Some general features of cosmic ray acceleration are summarized 

On ■ along with some inferences that can be drawn concerning the origin 

of the UHE component. The UHE luminosity density is found to be 

Qh| similar to that derived for GeV cosmic rays and its slope suggests a 

Q ■ distinct origin. Reports of clustering on small angular scale, if con- 

firmed, would rule out most proposed source models. More generally, 

j^ I it is argued that the highest energy particles can only be accelerated 

in sites that can induce an EMF £" ^ 3 x 10^'^ V and an associated 
power Lmin ^ £'^/Z ~ 10^^ W, where Z is the characteristic, electrical 
impedance, typically ^ lOOQ. Shock acceleration, unipolar induction 

C^ I and magnetic flares are the three most potent, observed, accelera- 

tion mechanisms and radio jet termination shocks, 7-ray blast waves, 
dormant black holes in galactic nuclei and magnetars are the least 
implausible, "conventional" manifestations of these mechanisms that 
have been invoked to explain the UHE cosmic rays. Each of these 
models presents problems and deciding between these and "exotic" 
origins for UHE cosmic rays, including those involving new particles 
or defects will require improved statistical information on the ener- 
gies, arrival times and directions, as should be provided by the Auger 
project. 



1 INTRODUCTION 

I have been asked to summarize "conventional" schemes for the acceleration 
of UHE cosmic rays, though any physical process capable of endowing a sub- 
atomic particle with the kinetic energy of a well-hit baseball/cricketball can 
hardly be considered conventional. This means that I shall leave others to 
review mechanisms that attribute the origin of these particles to topological 
defects, strings, monopole decay, supersymmetric hadrons, cosmic necklaces, 
cryptons and so on. Indeed, I suspect that the "hidden agenda" is for me 
to fail at my appointed task, and, like my colleagues on the MACHO exper- 
iment, to make the world safe for elementary particle theorists. I shall not 
disappoint. 

Many of the issues that I will cover have been recognised for some time 
and have been well-discussed in several excellent reviews including Hillas 
(1984) and Cronin (1996) and the many relevant contributions to the re- 
cent conference on this subject (Krizmanic, Ormes & Streitmatter 1998), 
including, especially, the lively summary by the late David Schramm. The 
conference proceedings edited by Chupp & Benz (1994) is also relevant. 

2 THE COSMIC RAY SPECTRUM 

In order to give this topic some context, consider the complete cosmic ray 
spectrum {eg Berezinski et al 1990). This extends over nearly twelve decades 
of energy from the proton rest mass, ~ 1 GeV, where their energy density 
is that of the microwave background, to at least 300 EeV (= 50 J = 3 x 
10~^ mpi). We can consider the cosmic ray spectral energy density inferred 
at the solar system, U{E) = {ATr/c){dI /d\nE) - 5 x 10-i^(E/10GeV)-°-^ J 
m~^ (correcting for solar modulation) extending from ~ 10 GeV to the "knee" 
at ~ 100 TeV- 10 PeV. (10 GeV cosmic rays are about 10 m apart and have 
an energy density comparable with that of the microwave background. The 
spectrum steepens above the knee: U{E) ~ 4 x 10^"'^^(ii^/10PeV)^^'"'^ J m^^. 
It then dips and flattens around the "ankle" (~ 1 — 10 EeV). UHE cosmic rays 
- the toenail clippings of the universe - are observed up to 300 EeV and, with 
a little imagination, U{E) ~ 1.5 x 10~^^(i?/10EeV)"*'-^ J m~^, comparable 
with the estimated, integrated background from 7-ray bursts. (Despite the 
large uncertainty, and the fact that the number density has fallen by ~ 10 



orders of magnitude, we do measure the EeV spectrum better than the MeV 
spectrum, of which we are, quite decently, ignorant.) 

The ~ 1 GeV-100 TeV cosmic rays are of Galactic origin. The ratio of Li, 
Be, B secondaries to C, N, O primaries measures their range to be X{E) ~ 
100(-E/10GeV)~'^'^kg m~^ {eg Axiord 1994). The cosmic ray luminosity of the 
Galaxy is then estimated as ~ MdU{E)c/X{E) ~ 2 x 10^^{E/lGeV)-^-^ W 
where M^ is the gas mass of the disk. Scaling from the local galaxy luminosity 
density (per Inii^ and assuming h = 0.6), we derive an average, cosmological, 
luminosity density (per \nE), < C> (E) r^Ax 10-''^^(^/10Gev)-°-^ W m'^, 
for lOGeV < E < lOOTeV. (For comparison, the stellar luminosity density 
is ~ 10^33 W m-3.) 

The UHE particles are almost surely extragalactic. As with 7-ray bursts, 
there is no good evidence for disk, halo, cluster or supercluster anisotropy 
(despite some tantalising hints in the past), (Takeda et al 1999). Further- 
more, magnetic confinement by the Galaxy is impossible - the Larmor radius 
rL{E) of a 300 EeV cosmic ray in a /iG field is ~ 300kpc. If we assume 
that UHE cosmic rays are protons, (and assuming that they are not, only 
makes matters worse), then they have a short lifetime to photo-pion produc- 
tion on the microwave background, (Greisen 1966, Zatsepin & Kuzmin 1966). 
The characteristic lifetime of a ~ 60 — 300 EeV cosmic ray is, very roughly, 
T{E) ~ 0.1(-E/300EeV)~^ Gyr. This implies that the luminosity density 
increases with energy < £ > (E) ~ U{E)/T{E) ~ 10-^'^{E/300EeVy-^ W 
m~^. At the highest measured energy, the estimated cosmological luminosity 
density is not significantly different from that of 10 GeV cosmic rays. The 
change in slope in the source spectrum, above ~ 60 EeV, is a strong indica- 
tion that these UHE cosmic rays comprise a quite distinct component from 
their lower energy counterparts. 

In order to investigate this further, it is necessary to take account of the 
fiuctuations in energy loss. Taking the 17 events reported by the AGASA 
collaboration above 60 EeV, it is possible to derive a maximum likelihood 
estimate of the unnormalized energy density, uncorrected for biases in de- 
tection efficiency. I find that if U{E) oc E°;E > Emm = 60 EeV, then 
a = —1.2 ± 0.5. I then calculate the probability that a particle of energy 
Eq has energy > E after time t, P{E,t; Eq), following Aharonian & Cronin 
(1994) and Bahcall & Waxman (1999). If we assume a power law for the 
luminosity density jC{E) oc E^; E > -Emin = 60 EeV, then the logarithm of 



the likelihood for obtaining the observed events is 
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Maximizing this function with respect to variation of (3, gives the estimate 
f3 = 0.3 ± 0.2. A more sophisticated computation that takes into account 
the detection probabilities of the different events should be performed, but 
it is unlikely to change the conclusion that the spectral luminosity density 
actually increases in the 60-300 EeV energy range and may even be consistent 
with a single, "top down" source with energy well above 300 EeV. 

There have been reports that UHE cosmic rays are significantly clustered 
on the sky. Specifically, in a sample of 47 events observed with AGASA, 
there are three pairs and one triple above 40 EeV with separations < 2.5°, 
comparable with the positional errors (Takeda et al 1999). (There are two 
more coincidences with events drawn from other samples.) There is no clear 
pattern for the associated particles to be ordered in energy and, in particular, 
one double has a 106EeV particle arriving over 3 yr. after a 44 EeV particle. 

// these associations are real, then there are three important implica- 
tions. Firstly, as particles are likely to be defiected by intergalactic mag- 
netic field through an angle 56 ~ {D^bY /^biE) then they will be delayed 
by ~ D"^ a. B / T l{EY c oc -E~^, where £b is the field correlation length and 
D ~ cT(150EeV) ~ 30 Mpc is the supposed source distance ( cf Miralda- 
Escude & Waxman 1996). Even Aesop would be challenged to explain how a 
~ 40 EeV cosmic ray precedes a ~ 100 EeV cosmic ray if they started at the 
same time and we must conclude that the source persists for several years, 
at least. This would rule out all particle/defect and 7-ray burst models. Sec- 
ondly, the small defiection angles at low energy limit the intergalactic field 
strength to -B < 20(£b/1Mpc)~^/^ fT, far smaller than generally supposed, 
though probably not excludable by direct observation. Thirdly, the presence 
of three ~ 40 EeV cosmic rays associated with high energy cosmic rays of 
much shorter range, implies that the background of low energy cosmic rays 
not associated with high energy events must be larger than the incidence of 
clustered events by roughly the ratio of their typical lifetimes ~ 30, which 
more than accounts for the remainder of the low energy sample. This, in 
turn, implies that the high energy cosmic rays must come from a very few 
sources which are, consequently, quite energetic: E ~ 10''^(r/3yr) J, where r 



is their lifetime. (If the low energy cosmic rays are scattered through ~ 2.5°, 
then r ^ 10^ yr and E ~ 3 x 10^^ J.) 

However, this clustering hypothesis, which is necessarily a posteriori when 
expressed in detail, is only supported with modest confidence. (A simple, 
Monte Carlo simulation distributing 47 points at random on half the sky 
and looking for similar patterns is quite instructive.) Particle/defect/burst 
explanations of UHE cosmic rays need not yet be rejected on these grounds. 

3 COSMIC RAY ACCELERATION 

The standard model of bulk cosmic ray production is first order Fermi accel- 
eration at strong, super- Alfvenic, shocks associated with supernova remnants 
and, possibly, winds from hot stars {eg Blandford & Eichler 1987). A typical 
relativistic proton will cross a shock, travelling with speed m, 0{c/u) times, 
gaining energy AE/E = 0{u/c) each traversal through scattering by hydro- 
magnetic waves moving slowly with respect to the converging fluid flows on 
either side of the front. The net mean relative energy gain is 0(1), but the 
process is statistical and a kinetic calculation shows that the transmitted 
spectrum will be a power law in momentum, f{p) oc p~'^^l^'^~^\ where r, (= 4 
for a strong shock), is the compression ratio. This mechanism can account, 
broadly, for the power (eg Malkov 1999), the slope {eg Axford 1994) and the 
composition {eg Ellison et al (1997)) of GeV cosmic rays. Shock acceleration 
is also, arguably, observed directly in SN1006 (Koyama et al 1995), as well 
as in the solar system {eg Erdos & Balogh 1994). 

The maximum energy to which a particle can be accelerated at a shock 
front is dictated by the scattering mean free path, (i{E) ~ {B/SB)l^rL{E), 
where SB is the amplitude of resonant hydromagnetic waves with wavelength 
matched to the particle Larmor radius. The diffusion scale-length of cosmic 
rays ahead of the shock is ~ £c/u and, assuming that this is limited by 
the size of the shock ~ i? we arrive at the unsurprising result that the 
maximum energy achievable in shock acceleration, assuming 6B < B and 
the presence of a large scale magnetic field, is -Emax = e£ ^ euBR ~ ed^/dt, 
the product of the charge and the motional potential difference across the 
whole shock. Equivalently, we conclude that in order to accelerate a proton 
by this mechanism to an energy £e, the rate of dissipation of energy exceeds 
-^min ~ £"^12, where Z = iiqu = iiqE / B is the effective impedance of the 



accelerator in SI units. 

Imposing this condition for a supernova remnant in the interstellar medium 
leads to an estimate -Emaxsnr ~ 30 TeV {eg Axford 1994), close to the knee. 
An additional source is needed between the knee and the ankle, where the 
source is generally supposed to be metagalactic. Larger shocks, especially 
those at Galactic wind termination shocks {eg Jokipii & Morfill 1987) and 
associated with gas flows around groups and clusters of galaxies have been 
invoked {eg Norman et al 1995). These shocks are likely to be relatively 
weak and therefore to transmit steeper spectra, as observed. The major un- 
certainty is the strength of the magnetic field. If i? ~ 30 pT at a galactic 
shock and ~ 10 pT at a cluster shock, then Emax ~ 10, 100 PeV respectively. 
Neither site is likely to accelerate the highest energy particles. 

An alternative accelerator is the unipolar inductor {eg Goldreich & Ju- 
lian 1967). The archetypical example is a pulsar - a spinning, magnetised, 
neutron star. The surface field will be quite complex but a certain quantity 
of magnetic flux <l> can be regarded as "open" and tracable to large distances 
from the star, (well beyond the light cylinder). As the star is an excellent 
conductor, an EMF will be electromagnetically induced across these open 
field lines £ ~ 0$, where $ is the total, open magnetic flux. This EMF 
will cause currents to flow along the field and as the inertia of the plasma 
is likely to be insignificant the only appreciable impedance in the circuit is 
related to the electromagnetic impedance of free space Z ~ 0.3/ioc ~ 100 Vt. 
The maximum energy to which a particle can be accelerated is -Emax ~ e£ 
and the total rate at which energy is extracted from the spin of the pulsar 
is Lmin ~ S"^ / Z. Taking the Crab pulsar as an example, -Emax ~ 30 PeV for 
protons and Lmin ~ 10^^ W. As the stellar surface may well comprise iron, 
even the Crab pulsar has the capacity to accelerate up to EeV cosmic rays. 
However, it is not obvious that all of this potential difference will actually 
be made available for particle acceleration. In particular, this is unlikely to 
happen in the pulsar magnetosphere as a large electric field parallel to the 
magnetic field will be shorted out by electron-positron pairs, which are very 
easy to produce, and radiative drag is likely to be severe. A more reasonable 
site is the electromagnetic pulsar wind and the surrounding nebula where 
particles can gain energy as they undergo gradient drift between the pole 
and the equator (Bell 1992). Pulsars may well contribute to the spectrum of 
intermediate energy cosmic rays. 

A third, protoypical accelerator is a fiare, for example one occuring on the 
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solar surface or the Earth's magnetotaiL Here magnetic instabihty leads to a 
catastrophic rearrangement, which must be accompanied by a large inductive 
EMF. Unless the instabilities are explosive, the effective impedance is again 
~ fiQU, where -u is a characteristic speed. Non-relativistic flares generally 
convert most of the dissipated magnetic energy into heat and are notoriously 
inefficient in accelerating high energy particles. 

Other acceleration mechanisms have been proposed and may contribute to 
the acceleration of the bulk of Galactic cosmic rays and relativistic electrons 
in non-thermal sources. These include a variety of second order processes and 
steady, magnetic reconnection. Many of them can be observed to operate 
within the solar system. However, they are thought to be too slow to be 
relevant to the acceleration of the highest energy cosmic rays. 

4 ZEVATRONS 

Having argued that the three most potent, observed acclerators are shocks, 
unipolar inductors and flares, let us see how they can be modified to account 
for ~ ZeV cosmic rays. Firstly, note that, as m ~ c, mildly relativistic 
shocks minimise the power that has to be invoked to attain high energy. 
Specifically, we need a power > 10^^ W to account for 300 EeV cosmic rays 
and this exceeds the bolometric luminosity of a powerful quasar. One of 
the few sites where such a large potential difference can be achieved is the 
termination shock of a powerful radio jet like that associated with Cygnus 
A {eg Cavallo 1978). Stretching the numbers a little, we combine a field 
strength ~ 10 nT, with a speed ~ c and a transverse scale ~ 3 kpc which 
gives -Emax ~ 300 EeV. The problem with this model is that observed UHE 
cosmic rays are not positionally identified with the few known radio sources 
within D{E) ~ 30 Mpc that might be powerful enough to account for them 
(c/Farrar & Biermann 1998). 

A more elaborate shock accelerator is the 7-ray burst blast wave {eg Wax- 
man 1995). Here, the shocks (assumed to be spherical) are ultrarelativistic 
with Lorentz factor F. The maximum energy, measured in the frame of the 
explosion, to which a proton can be accelerated in a dynamical timescale 
from an ultrarelativistic shock of radius R is -Emax ~ eB'Rc, where B' is the 
comoving field strength. The explosion power, adopting the most elemen- 
tary of assumptions, is then Lmin ~ 47rF^(_Einax/e)^//xoc. Observed bursts 



have typical explosion powers estimated to be Lcxp ~ 10^^ W, which can be 
consistent with 300 EeV proton acceleration as long as F < 300, which is 
just compatible with existing models. A serious physical constraint is the 
avoidance of radiative loss in this environment. An observational concern 
with this model is the improbability of having enough active bursts close 
to supply the highest energy particles roughly isotropically (c/Waxman & 
Miralda-Escude 1996). 

The most relevant variant on unipolar induction is magnetic energy ex- 
traction from spinning, black holes, where the magnetic field is supported 
by external current, and the horizon is an imperfect conductor with resis- 
tance ~ 100^2 (eg Thorne it et al 1986). This impedance is matched to the 
electromagnetic load so that roughly half of the available spin energy ends 
up in the irreducible mass of the hole, the remainder being made available 
for particle acceleration. The total electromagnetic power needed to account 
for ~ 300 EeV acceleration is, once more, ~ 10^^ W. A rapidly spinning, 
~ 10^ Mq hole endowed with a field strength ^ 1 T or a ~ 10^ M© hole 
threaded by a ^ lO'* T field suffices to accelerate 300EeV particles. The 
major concern with this model is that the radiation background must be ex- 
tremely low in order that catastrophic loss due to pion and pair production 
be avoided. Specifically, it is necessary that the microwave luminosity in an 
acceleration zone, of size R, be < 10^'*(-R/10^^m) W, far smaller than the 
unobserved electromagnetic power. 

The best generalization of fiare acceleration involves "magnetars" which 
are young, spinning neutron stars endowed with a ~ 10 — lOOGT surface 
magnetic field as first postulated by Thompson & Duncan, (1996). Now, 
the observation of 5 — 7 s period pulsations from three "soft gamma re- 
peaters" effectively confirms their identification as old magnetars that have 
been decelerated by electromagnetic torque and which are now powered by 
magnetic energy which is released in a series of giant fiares, Kouveliotou et 
al 1998. The inductive EMFs associated with an electromagnetic fiare from 
a magnetar can be as high as £ ~ 3 x 10^^ V, making them candidate UHE 
accelerators because the surface composition is likely to be Fe. However, the 
available reservoir of magnetic energy is only ~ 10^° J and the magnetar 
birthrate is no more than ~ 10^^ yr^^ in the Galaxy. This rules them out 
as an extragalactic source. Only if UHE cosmic rays have a Galactic origin, 
(and the large scale anisotropy observations suggest quite strongly that do 
not), can there be enough power in magnetars to account for the UHE energy 



density. 

5 DISCUSSION 

I have argued, tentatively, that UHE cosmic rays are created in a new pop- 
ulation of extragalactic sources with an average luminosity density that ap- 
proaches that of Galactic cosmic rays. I have also described problems with 
each of the candidate "conventional" mechanisms for accelerating protons 
to these high energies. Quite different, general inferences have been drawn 
here, from the same data, by Waxman, and elsewhere by others. All of this 
underscores the need for better statistics which should be met by the Auger 
project. Perhaps the most pressing need is to understand if particles of very 
different energy have a common origin. If true, this must rule out essentially 
all primordial particle/topological defect, neutron star, 7-ray burst explana- 
tions, leaving only massive black holes and radio source models among the 
possibilities discussed above. In this case, it will be possible to seek identifi- 
cations, especially at the highest energies, where the positions will be most 
accurate and the delays due to magnetic scattering the smallest. If, alter- 
natively, clustering and its implications are not substantiated, then the next 
best clues will probably come from composition studies and detailing the 
large scale distribution on the sky. 

The most exciting outcome of all of this is that we are dealing with a new 
particle or defect with energy well out of the range of terrestrial accelerators. 
(For example, if there is a particle of energy Ex which decays with half life 
Tx into A^ protons, then the cosmological energy density of these particles 
must be fix ~ 3 X 10-^N-\Ex/lYeV)-\HoTxy\) Whatever happens, 
in a subject where the dullest and most conventional theories involve mas- 
sive, spinning, black holes, ultrarelativistic blast waves and 100 GT fields 
threading nuclear matter, the future is guaranteed to be interesting. 
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